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Abstract 

The Search for Extraterrestrial Intelligence (SETI) has thus far failed to provide a convincing detection 
of intelligent life. In the wake of this null signal, many "contact pessimistic" hypotheses have been for- 
mulated, the most famous of which is the Rare Earth Hypothesis. It postulates that although terrestrial 
planets may be common, the exact environmental conditions that Earth enjoys are rare, perhaps unique. 
As a result, simple microbial life may be common, but complex metazoans (and hence intelligence) will 
be rare. This paper uses Monte Carlo Realisation Techniques to investigate the Rare Earth Hypothesis, 
in particular the environmental criteria considered imperative to the existence of intelligence on Earth. 

By comparing with a less restrictive, more optimistic hypothesis, the data indicates that if the Rare Earth 
hypothesis is correct, intelligent civilisation will indeed be relatively rare. Studying the separations of 
pairs of civilisations shows that most intelligent civilisation pairs (ICPs) are unconnected: that is, they 
will not be able to exchange signals at lightspeed in the limited time that both are extant. However, the 
few ICPs that are connected are strongly connected, being able to participate in numerous exchanges of 
signals. This may provide encouragement for SETI researchers: although the Rare Earth Hypothesis is in 
general a contact-pessimistic hypothesis, it may be a "soft" or "exclusive" hypothesis, i.e. it may contain 
facets that are latently contact-optimistic. 
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1 Introduction 

The attributes of the planet Earth are of critical importance to the existence and survival of life upon it. In 
fact, it may be so finely tuned that few planets in the Galaxy share its life-friendly characteristics. From this 
premise, it is almost inevitable to reach the conclusion that intelligent life (at least any that is predicated on 
evolution from complex metazoans) is also rare - perhaps unique to the planet Earth. 



These ideas have been encapsulated in what is known as the Rare Earth Hypothesis (Wa rd & Brownlee . 
2000). It can be summarised thus: 

1 . Simple life may be commonplace in the Universe. The existence of extremophilic organisms in what 
were originally considered to be inhospitable regions (hydrother mal vents, acidic pools, toxic waste, 



deep in the Earth's crust) has shown the hardiness of simple life dCavicchiol"ill2002l : lDiaz & Schulze-Makuch 



2006). Inde ed, these habitats are believed to be dup licated els ewhere in the Sol ar Syste m - e.g. Mars 
dFormisano et all 12004 iKrasnopolskv et all |2004. Europa dCarr et all 1 19981) . Titan dStofan etal. , 
20071) . Enceladus dParkinson et all l2007t llSpencer & Grinspooni 120071) - so it is still possible that 
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"alien" life may be closer to home than once thought. 

However, although simple life is resilient and adaptable, the evolution of complex animal life is 
extremely difficult. For this to be achieved, there are certain criteria (hereafter referred to as the 
Earth Criteria) that must be satisfied, in order for animals to thrive. 



A (non-exhaustive) list of the Earth Criteria follows: 



A pla net within a critical range of orbital radii - the "stellar habitable zone" (Hart, 1979; K asting et al 



1993) 



• A star within a critical mass range (large enough to push the habitable zone outside the planet tidal 
locking radius, and small enough to provide sufficient energy while avoiding UV exposure) 



• A sta r located in a critical region of the Galaxy - the "galactic habitable zone" dLineweaver et al 
2004 



• A planet within a critical mass range to maintain a suitable atmosphere 

• A planet with a stable low eccentricity orbit (to avoid extreme tempera ture changes). This also 
requires a relatively large moon to provide axial stability d WalthamL 12004 ) . 

• A planet with sufficient raw materials to generate amino acids and proteins 

• A planet with suitable atmospheric composition, in particular the production of atmospheric oxygen 
- initially produced by cyanobacteria in Earth's early history (ICanfieldl 120051) 



A planet with plate tectonic activi ty to regulate atmospheric composition and the balance of carbon 
(IBounamam Christine et all [2007)) 



• The pres ence of Jupiter to control the rate of comet and a steroid impacts - although this is now in 
question dHorner & Jonesl l2008i 120091: iHorner et alll2010l) . 



The weakness of this hypothesis rests in the (usually implicit) assumption that all the Earth Criteria are inde- 
pendent of each other. Taking Jupiter as an example: asking whether Jupiter exists or otherwise in the Solar 
System is not meaningful, as planet formation is a complex, non-linear process: every planet in the Solar 
System owes its formation to its surounding environment, and therefore its planetary neighbours, through 
the d ynamics of migration dRavmond et all 120061: Paardekooper & P apaloizou. 2008) planet-planet scatter- 
ing dFord & RasioL 120081: iRavmond et alll2009h . resonances dCresswell & Nelsonll2006l) . and other secular 
phenomena (e.g. iBatvgin & Laughlinl l2008). Without Jupiter, the Earth as it is today may not have formed 
at all. 



Testing the Rare Earth Hypothesis 



3 



This paper will investigate the influence of a subset of the Ea rth Criteria on the re s ulting distribution of 
inhab ited planets, using Monte Carlo Realisation Techniques dVukotic & Cirkovic , 2007 , 2008 : iForganl 
2009). The paper is structured as follows: section [2] will outline the methods used to simulate the distribu- 
tion of life in the Galaxy; section[3]will define the input parameters for the simulations run; the results are 
displayed in section|4]and summarised in section|5] 



2 Method 



The n umerical simulations are carried out using the Monte Carlo Realisation techniques outlined in Forganl 
(l2009h . hereafter referred to as Paper I. A brief summary of the method follows for completeness. 

In essence, the method generates a Galaxy of N* stars, each with their own stellar properties (mass, lu- 
minosity, location in the Galaxy, etc.) randomly selected from observed statistical distributions. Planetary 
systems are then generated for these stars in a similar manner, and life is allowed to evolve in these planets 
according to some hypothesis of origin. The end result is a mock Galaxy which is statistically represen- 
tative of the Milky Way. To quantify random sampling errors, this process is repeated many times: this 
allows an estimation of the sample mean and sample standard deviation of the output variables obtained. 

The inputs used to define the mock Galaxy (e.g. the Galaxy's surface density profile, the initial stellar mass 
function (IMF), the star formation history (SFH), etc.) are of critical importance. Paper I focussed on using 
current empirical data (especially for the simulation of exoplanets) to define the mock Galaxy. This paper 
will attempt to improve on the inputs of Paper I: these improvements are discussed below. 

2.1 Improvements on the Model 
2.1.1 The Simulation of the Galaxy 

In Paper I, the Milky Way was simulated in two dimensions only (in polar coordinates (r, </>)). As a first 
improvement, the Gala xy is given vertical structure, incorporating both the thick and thin stellar discs 
dOsthe&CarrollLll996l) : 



p(r, z) = n e- r ^' /rH f e -W**Mn + QMe-^'^A (1) 
Secondly, the metallicity gradient of the Milky Way was previously simulated using only one curve: 

Z, = -z grad \ogl^-\ (2) 
In tru th, there are many differing measurements of the abundance gradient in the Galaxy dRolleston et al 



2000), dependent on the metals studied. This reflects the different synthesis processes at work for differing 
elements. This can be (crudely) reproduced by allowing z gra d to have a distribution of values - in this case, 
a Gaussian distributio n, with sample mean and sample standard deviation defined by the measurements of 
iRolleston et~ai]d2000l) . 

Finally, measures h ave been taken to c orrela te the age and metallicity of the stars. The Age Metallicity 
Relation (AMR) of iRocha-Pinto et al.l d2000l) (with its errors) defines upper and lower bounds to the age 
of a star (given its metallicity). The star formation history (SFH) has also been improved, allowing better 
time resolution ( Rocha-Pint o et al . 2000, see Figure Q]). 



2.1.2 The Simulation of Stars 

An important change to this work is the sim ulation of stellar luminosity evolu tion. As stars evolve along the 



main sequence, their luminosity increases (Schro der & Connon S mith. 2008J). As the luminosity increases 



the lo cation of the stellar habitable zone must move further away from the star (Hart, Il979t iKasting et al 



1993). This implies that any planets with biospheres initially in the habitable zone can leave the habitable 
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Age (Gyr) 



Figure 1 : The star formation history used in this work nRocha-Pinto et ai, 2000) 



zone on a timescale thz> which is a function of the host star's initial luminosity and the planet's distance 
from it. Together with the main sequence lifetime Tms> it defines a maximum lifetime for any biosphere: 



MIN{tms,t HZ ) 



Where tms is given by 



tMS 
*MS,0 



1 -3 



Mr. 



(3) 



(4) 



For more information, see Prialnik (2000) (Note that quantities referring to the Sun have the subscript 0). 
T he luminosity evolution of the stars are approximated by extrapolating the simulated solar luminosity data 
of ISchroder & Connon Sm ith (2008) to all main sequence starsj: 



L(t) = 0.7 + 0.144 



Gyr 



(5) 



2.1.3 The Simulation of Planets 

Current exoplanet data, while impres sive, is still incomplete. This introduces significant bias into the re- 
sults of any simulation dForganl 120091) . and precludes the discussion of the Rare Earth Hypothesis if using 
observations alone (as statistical analyses cannot currently simulate Earths with any robustness). To bypass 
this problem, the empirical data is replaced by theoretical relations: this allows the simulation of planetary 
objects down to lunar masses. 

The probabilit y of a star hostin g planets is a function of its metallicity: this code uses the distribution as 
described bv lWvatt et ail d2007h : 



P(z) = 0.03 x 10 z s 
The Planetary Initial Mass Function (PIMF) is approximated by a simple power law 

(Mp)- 1 



(6) 



(7) 



P{M P 

which operates over the mass range of [M moon , 25 M j up ] . To correctly reproduce the distributio n of plane- 
tary radii, two different radii distribution functions are used. Jovian planets reproduce the data of Armitape 



1 This may seem a weak assumption, but the stars of interest in these simulations will be close to solar type, so the approximation 
is reasonable in this first instance. 
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Figure 2: The mass-radius relation for a sample of planets. 



(200 7]), wh ich accounts for the effects of Type II planetary migration. For terrestrial planets, the data of 
Ida & Linl d2008l) (Fig. lc) is emulated: a simple parametrisation allows the trend for low mass objects to 
be recovered (see Figure O. It should be noted that in essence this is swapping one weakness for another: 
while the bias of empirical data is lost, the uncertainty of current planet formation models is gained. 



Also, as the mass function can simulate Moon-mass objects, any object with a mass less than Pluto's 
that resides within another planet's Hill sphere (that is, it resides within the gravitational influence of said 
planet) is considered a moon of that planejl This property will become important for the simulation of the 
Rare Earth Hypothesis. 



2.1.4 The Simulation of Life 

The simulation of life proceeds in much the same manner as Paper I: life must achieve several difficult 
goals, e ach goa l req uiring a time r, to be achieved, in order to evolve into a technologically capable, sentient 
species (ICarten. 120081). During this process, resetting events can occur (e.g. cometary impacts, supernovae, 
gamma ray bursts, see lAnnisll999l) . which cause large scale extinctions, and reset the evolution of life to an 
earlier stage (or in the worst case, sterilise the planet entirely). The reader is referred to Paper I for a more 
in depth discussion of how this is achieved in the code. A significant change to the model is the mechanism 
for resetting events becoming sterilisation events. This was described by a probabi lity of annihilatio n 



Pannihiiate, which was created in order to define the Galactic Habitable Zone dLineweaver et all 12004). 



This simplified approach has been replaced with an attempt to simulate the loss of biodiversity incurred 
as the result of a reset. The effect of a reset is to reduce biodiversity by a fraction x: if x is greater 
than 1, the planet is sterilised. Appealing to the Central Limit Theorem, x is selected from a Gaussian 
Distribution with mean 0.5, and standard deviation 0.25. This means that on average, 5% of resets will 
result in annihilation. The average number of resets increases with proximity to the Galactic Centre - this 
is now parametrised by: 



presets — l^Earth^ 



\ r gal,@ J 



(8) 



Where HEarth = 5 dRaup & Sepkoskil 119821) . This provides the mean of a Gaussian distribution from 
which N resets is sampled. The habitation index of Paper I (which is assigned to all planets in the simula- 
tion) is modified to account for the potential existence of microbial life: 



2 This does not guarantee the observational rule of thumb that the probability of a terrestrial planet having a substantial moon is 
around 0. 25. However, this is a reasonable first approximatio n, with more detailed studies of this issue requiring future observational 
input (e.g. lSartoretti & Sc hneider 1999; Kipping et alj2009l) . 
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Biosphere which has been annihilated 

Planet is lifeless 

Planet has microbial life 

Planet has primitive animal life 

Planet has intelligent life 

Planet had intelligent life, but it destroyed itself 

Planet has an advanced civilisation 



(9) 



Note that interplanetary colonisation is not modelled in this work. 

2.2 New Outputs - Civilisation Interaction 

As the code produces data pertaining to each individual civilisation, it is possible to study the entire dataset 
for each run, and identify the potential for communication between all possible pairs of civilisations. For 
N galactic civilisations, there are N(N — l)/2 pairs of civilisations. For each intelligent civilisation pair 
(ICP), the following outputs can be calculated: 

1 . Their physical separation in kiloparsecs (dx), 

2. The available window of communication dt (that is, the maximum time interval where both civilisa- 
tions exist and are able to communicate), 

3. The space time interval ds 2 — c 2 dt 2 — dx 2 . This quantity determines whether a signal travelling 
at lightspeed can traverse the distance between two civilisations within the communication window 
(assuming the intervening space to be Minkowskian). If ds 2 < 0, then the signal will fail to reach its 
destination before the window closes. If ds 2 = 0, then the signal will reach its destination at the same 
instant the window closes. If ds 2 > 0, then the signal will reach its destination within the window, 
and it is therefore possible for communication between the two civilisations to be established. 

4. The "contact factor" f contact = which counts how many "conversations" (pairs of signals) can 
travel between the two civilisations. 

These outputs give extra information on the distribution of civilisations in the Galaxy, and their potential 
connectedness by signals travelling at lightspeed. 

3 Inputs 

Two separate hypotheses were tested with this model. Each was subjected to 30 Monte Carlo Realisations 
(MCRs), with each realisation containing N star s = 10 9 . This is of course two orders of magnitude short 
of the Milky Way's stellar content, but computational constraints prevented increasing N stars further. The 
interested will be able to multiply subsequent results by 100 to obtain an estimate of Milky Way figures. 
In any case, absolute numbers are less relevant to the issue at hand: this study focuses on comparing two 
hypotheses, and comparing relative trends (which is a more reliable route in studies of this nature). 

3.1 The Baseline Hypothesis 

This basic hypothesis requires only that a planet must be in the stellar habitable zone for life to form upon it. 
If the planet's surface temperature lies between [0, 100]°C, then microbial life can for m upon it. Complex 
anima l life will only form if the planet's surface temperature lies between [4, 50]°C dWard & Brownlee, 
20001) . This hypothesis was tested to provide a comparison with the results of the Rare Earth Hypothesis. 



Testing the Rare Earth Hypothesis 



7 



10'° - 
10 s - 
10 6 - 
10 4 - 
10 2 - 
10°_ 

-2 4 6 - 4 

Habitation Index Habitation Index 

Figure 3: The habitation index for the Baseline Hypothesis (left) and the Rare Earth Hypothesis ( right). 

3.2 The Rare Earth Hypothesis 

This hypothesis builds on the baseline hypothesis by also requiring that animal life will only form on a 
planet if the following four conditions are met: 

1. The planet's mass is between [0.5, 2.0]M^, 

2. The star's mass is between [0.5, 1.5]M©, 

3. The planet has at least one moon, (for axial stability and tides) 

4. The star system has at least one planet with mass > 10M ffi in an outer orbit (for shepherding aster- 
oids) 

4 Results & Discussion 

4.1 The Distributions of Life and Intelligence 

The properties of the Rare Earth hypothesis galaxy can now be compared against those of the Baseline 
hypothesis. Comparing the habitation index for both hypotheses (Figure |3), it can be seen that (by con- 
struction) microbial life (linhabit = 0.5) is unaffected by the Rare Earth hypothesis, whereas the prevalence 
of animal life (linhabit — 1) is reduced by a factor of 10 4 against the baseline. This reduction is thus prop- 
agated into the intelligent biospheres (linhabit > 2). However, despite some quite stringent conditions 
on the planetary system architecture (conditions 3 and 4, section l3~2l the number of intelligent biospheres 
numbers in the hundreds: the implications for SETI are discussed in the next section. 

As stellar mass is a key condition to the Rare Earth Hypothesis, it should be expected that the two hypothe- 
ses' distributions diverge, and this is indeed the case: Figure 0] shows the distribution of stellar mass for 
both hypotheses. The IMF is modified by the effects of the habitable zone (and the distribution of exoplanet 
semi-major axis) to give the characteristic bump between 1 and 2 solar masses. Comparing the hypotheses 
shows that although the Baseline hypothesis favours lower mass stars for intelligent biospheres (for their 
increased longevity), the Rare Earth Hypothesis must discard the substantial number of stars that are less 
than 0.5M Q . 

This bias towards lower mass should be reflected in the distribution in semimajor axis (as lower mass stars 
have closer, more stationary habitable zones). Figure [5] shows that this is true for the Baseline Hypothesis 
(with intelligent biospheres dropping off as R > 1.5 AU), and doubly true for the Rare Earth Hypothesis, 
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Figure 4: Stellar Mass for the Baseline Hypothesis (left) and the Rare Earth Hypothesis (right). The black 
lines indicate all biospheres, the blue lines indicate all intelligent biospheres. 
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Figure 5: Planet semimajor axis for the Baseline Hypothesis (left) and the Rare Earth Hypothesis (right). 
The black lines indicate all biospheres, the blue lines indicate all intelligent biospheres. 



selecting a narrow radial range between [0.8, 1.9] AU. 

The most striking difference can be seen in the distribution of galactocentric radius (Figure [6]). While the 
Galactic Habitable Zone (GHZ) can be identified in the Baseline Hypothesis (with a small contingent at 
lower radii, which presumably exists due to the lack of modelling of the central supermassive black hole 
(SMBH) and hypervelocity stars in the inner regions), the Rare Earth Hypothesis appears to have no GHZ. 
This is unexpected: the four conditions of the Rare Earth Hypothesis tested here do not affect where intel- 
ligent systems should lie; why then does the GHZ not appear (with reduced numbers)? 



4.2 Communication and Connectivity 

The prima facie conclusion (having studied the results of the previous section), is that if the Rare Earth 
Hypothesis is correct, and intelligent civilisations are infrequent, then the potential for communication is 
also low. This expectation can be tested by calculating the interaction variables discussed previously. As 
the focus has now shifted from individual intelligent civilisations to intelligent civilisation pairs (ICPs), 
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Figure 6: Galactocentric radius for the Baseline Hypothesis (left) and the Rare Earth Hypothesis (right). 
The black lines indicate all biospheres, the blue lines indicate all intelligent biospheres. 
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Figure 7: Separations of ICPs for the Baseline Hypothesis (left) and the Rare Earth Hypothesis ( right). 



the numbers duly increase 

from N t0 M^-i2. Fi 

gure [7] shows the distribution of ICP separation dx for 
both hypotheses. The baseline hypothesis exhibits a sharp peak at around 8 kpc (the location of the GHZ), 
accompanied by a long decay. This distribution is re miniscent of the lognormal distribution expected if the 
tools of the Statistical Drake Equation were applied (Maccone, 20091) . The distribution reaches its mode in 
steps: these steps are presumably sensitive to the local galactic spiral structure. The Rare Earth Hypothesis 
has no apparent GHZ, so the distribution (though reduced in magnitude) peaks at a much lower 3 kpc. 



Does this reduced separation imply increased connectivity? The answer depends on the communication 
window for the ICP. The longer the window (i.e. the larger overlap in history where both civilisations 
exist), the longer the separation can be while allowing the ICP to be connected. The Baseline Hypothesis 
favours shorter communication windows (Figure [HJ, which reduces the connectivity. Apart from small 
fluctuations at larger values, the Rare Earth Hypothesis agrees. When considering the space-time interval 
ds 2 (Figure [9), the reduced connectivity becomes apparent. ICPs that are unconnected (negative values) 
are much more frequent than connected civilisations (positive or zero values). This does not spell the end 
for SETI, however, when the contact factor (i.e. the number of conversations) is considered: although few 
ICPs enjoy the privilege of contacting each other, those that do can expect a great deal of conversation 
(Figure [Tol). where each hypothesis agrees that a select few will enjoy potentially thousands of exchanges 
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Figure 8: Communications Window (maximum time interval for communication) for ICPsfor the Baseline 
Hypothesis (left) and the Rare Earth Hypothesis ( right). 
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Figure 9: Space-time interval (ds 2 ) of ICPs for the Baseline Hypothesis (left) and the Rare Earth Hypoth- 
esis (right). Unconnected ICPs have ds 2 < 0, connected civilisations have ds 2 > 0. 



with other civilisations. 



4.3 Future Improvements to the Model 

Numerical modelling of this type is generally a shadow of the entity it attempts to model, in this case the 
Milky Way and its constituent stars, planets and other objects. While a substantial improvement over the 
work of Paper I, there is still potential for future work. Several suggestions are listed here. 

1 . A more accurate Galactic model, taking into better account its chemical diversity, stellar clustering 
and the inner regions (specifically the central supermassive black hole and the hypervelocity stars 
orbiting it). 

2. An improved planetary architecture model, better equipped to deal with moons and the planet mass- 
semimajor axis distribution. Also missing is the modellin g of orbital eccentricity and inclination, 
potentially of great importance in issues of habitability (e.g. IWilliams & Pollardll2002l ISpiegel et al.1 



2008). 
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Figure 10: Contact Factor (number of exchanged signal pairs) for ICPs for the Baseline Hypothesis (left) 
and the Rare Earth Hypothesis ( right). 



3. Improved modelling of the connectivity of civilisations (potentially extending to the modelling of 
interstellar colonisation and face-to-face contact). 



5 Conclusions 



This paper has tested the Rare Earth Hyp othesis ( Ward & Brownlee , 2000h using the Monte Carlo Realisa- 
tion techniques outlined in Forgan ( 20091) . By comparing the results to a baseline hypothesis, the influences 
of the criteria for a planet to be officially designated as "an Earthlike planet" can be studied. In this work, 
the criteria were limited to planet mass, star mass, the presence of a moon, and the presence of a Jupiter 
type object in a more distant orbit. It is shown that these criteria alone greatly reduce the number of in- 
telligent civilisations in the Galaxy (compared to the baseline). As expected, the stellar mass criterion 
results in a narrow range of planet semi-major axes where intelligent biospheres exist. Interestingly, the 
Galactic Habitable Zone, apparent in the Baseline Hypothesis, was not visible in the Rare Earth Hypothesis. 



This result is important for civilisation connectivity: reducing the civilisation separation means that, for a 
given time interval of communication, civilisations under the Rare Earth Hypothesis are able to exchange 
more signals than civilisations under more contact-optimistic hypotheses. The implications for SETI are 
somewhat mixed: while Earth may be much more likely to be a disconnected than connected civilisation, 
if it is connected, it can expect substantial conversation from other civilisations (while the Sun remains 
in the Main Sequence). Therefore, the Rare Earth Hypothesis (i n the f ormul a tion descr i bed in this work) 
is a "soft" or "exclusive" hypothesis (using the nomenclature of iBrinl dl983l) : ICirkovid (12009b ). in that it 
is not a completely contact-pessimistic hypothesis, but one that is contact optimistic for a small subset of 
civilisations in the Galaxy. 
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